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The phenomenon of quantum antiresonance (QAR), i.e., exactly periodic recur- 
rences in quantum dynamics, is studied in a large class of nonintegrable systems, the 
modulated kicked rotors (MKRs). It is shown that asymptotic exponential localiza- 
tion generally occurs for r/ (a scaled K) in the infinitesimal vicinity of QAR points r]Q. 
The localization length is determined from the analytical properties of the kicking 
potential. This "QAR-localization" is associated in some cases with an integrable 
limit of the corresponding classical systems. The MKR dynamical problem is mapped 
into pseudorandom tight-binding models, exhibiting dynamical localization (DL). By 
considering exactly-solvable cases, numerical evidence is given that QAR-localization 
is an excellent approximation to DL sufficiently close to QAR. The transition from 
QAR-localization to DL in a semiclassical regime, as r/ is varied, is studied. It is 
shown that this transition takes place via a gradual reduction of the influence of the 
analyticity of the potential on the analyticity of the eigenstates, as the level of chaos 
is increased. 
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I. INTRODUCTION 

The study of "quantum chaos", i.e., understanding the fingerprints of classical chaos in 
quantum mechanics |l|J^, has led to the discovery of a variety of new quantum-dynamical 
phenomena. Several such phenomena occur in time-periodic systems described by the general 
Hamiltonian 

H = Ho + HJ{t), (1) 

where Hq is some time-independent Hamiltonian, Hi is a perturbation, and f{t) is periodic 
with period T, fit + T) = f(t). In many cases, f{t) is chosen, for simplicity, as a periodic 
delta function, f(t) = Ar(t) = J2'^-oo^i'^ " giving the well-known class of "kicked" 
systems. Representative models in this class are the kicked rotor [3,^^|6|,|7|,|8|,P JT0|JTT|JT^ , the 
kicked Harper model and the kicked harmonic oscillator 

The quantum dynamics of time-periodic systems (|I]) is governed by their quasienergy 
(QE) spectrum (i.e., the spectrum of the one-period evolution operator). Different prop- 
erties of the QE spectrum lead to quantum-dynamical phenomena having, in general, no 
classical analogue. A classic example is the quantum suppression of chaotic diffusion in 
the kicked rotor (KR) accompanied by quasiperiodic recurrences p. An important 
interpretation of this phenomenon has been given by showing first that, in the angular- 
momentum representation, the QE eigenstates of the KR satisfy the equation describing a 
ID tight-binding model with pseudorandom disorder. This pseudorandomness is generic, as 
it exists for almost all (irrational) values of a scaled (dimensionless) h, which we denote here 
by 1]. It was found [|15| that in several interesting cases pseudorandom tight-binding models 



exhibit localization properties similar to those of truly random ones (Anderson localization) 
||16|| . Assuming the general occurrence of this localization, it follows that, generically, the 
QE eigenstates are exponentially localized in angular momentum and the QE spectrum is 
pure point. This localization in pseudorandom tight-binding models, equivalent to quantum- 
dynamical systems with nonintegrable classical counterparts, is called dynamical localization 



(DL) |T^. The quantum suppression of diffusion in the KR is an immediate consequence 
of DL. Despite the fact that DL has no classical analogue, there exists a remarkable and 
simple relation between the classical chaotic-diffusion coefficient D in the KR and the asymp- 
totic DL length ^ in the semiclassical regime (sufficiently small rj): ^ ^ D/2 For 
nongeneric, rational values of 77, there occurs in all the kicked systems the phenomenon 
of quantum resonance 0], i.e., the quadratic increase of the energy expectation- value with 
time. This phenomenon is due to an absolutely continuous QE spectrum, exhibiting a band 
structure. 

In this paper, DL is approached in the light of a different kind of phenomenon for systems 
(HI): exactly periodic recurrences. This phenomenon is defined, in general, by 

= e-'^ , (2) 

where U is the one-period evolution operator for ([^), e~*^ is some constant phase factor 
(a c number), and p is the smallest positive integer for which (^ is satisfied. Thus pT is 
the recurrence period. As it will become apparent in this paper, the phenomenon (0) may 
occur, in general, only for very special values of 77, and it is thus nongeneric. In fact, for 
the general class of systems introduced in this paper, it occurs precisely at values of t] = rjo 
corresponding to quantum resonances in the kicked systems. At the same time, this 
phenomenon, manifesting itself in bounded, periodic variation of expectation values, is dia- 
metrically opposite to quantum resonance. We shall therefore refer to (H) as to the quantum 
antiresonance (QAR) phenomenon. 

While this phenomenon is nongeneric, we show in this paper that, for a large class 
of nonintegrable systems, it is generally accompanied by a very interesting effect: In the 
immediate vicinity of QAR (infinitesimal 1] — rjo), there takes place asymptotic exponential 
localization with a pure-point QE spectrum. The existence of this "QAR-localization" is rig- 
orously established in the framework of a self-consistent approach, which allows for an exact 
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determination of the asymptotic localization length ^o- For rjo = 0, the QAR-localization is 
associated with an integrable limit of the classical Hamiltonian. On the other hand, values 
of r/o 7^ usually correspond to the strong quantum regime of a nonintegrable Hamiltonian, 
exhibiting chaotic diffusion. 

As in the case of the KR and other systems |l6|,|7|j8| JlO|JTT|l , we show that our class of systems 
can be mapped into tight-binding models with pseudorandom disorder (for t] sufficiently ir- 
rational). DL is then expected to occur generically in our systems. For rj infinitesimally 
close to rjQ, the tight-binding models are not defined, so that QAR-localization cannot be 
viewed, strictly speaking, as a kind of DL. However, if rj is sufficiently irrational and close 
to riQ, one expects DL to take place and to exhibit approximately the same features as those 
of QAR-localization. In particular, the DL length ^ should be well approximated by the 
QAR-localization length ^o, and one expects that ^ — > as ?7 — ?7o. We provide strong 
numerical evidence that this is indeed the case. Since can be determined exactly, this 
seems to be the first case where a DL length ^ can be found with arbitrary accuracy in 
nonintegrable systems. 

We also study the dependence of ^ on 77 — ?7o for 77 not very close to 770. This allows one to 
understand the transition from QAR-localization to DL in regimes basically different in na- 
ture from QAR. Values of rj sufficiently far from t]q may correspond to semiclassical regimes 
of local or global chaos. As already mentioned, the semiclassical regime of global chaos in 
KR systems is characterized by the approximate relation ^ ^ D /2 0,0,^]^. Calculations of 
^ for the KR were performed by applying the method of minimal Lyapunov exponent 
p!8| to an equivalent pseudorandom tight-binding model. The method is based on a finite 
truncation of the (generally infinite) vector of hopping constants, so that only the first 
neighbors are kept. Then ^ = limAr^oo 1/7jv; where 77V is the minimal Lyapunov exponent 
of a 2A^-dimensional symplectic map associated with the truncated model. We show that 
the truncated model has, for all A^, a well-defined dynamical equivalent exhibiting QAR. 
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Then, by studying numerically the dependence of on both and 77, we show that the 
transition from QAR-localization to DL in a semiclassical regime ~ D /2) takes place via 
a gradual reduction of the influence of the analyticity of the system on the analyticity of the 
eigenstates, as the level of chaos is increased. 



In the simple case of = 1, i.e., a nearest-neighbor pseudorandom "Lloyd model", we 
derive the exact relation ^0 = -D/2. Since many numerical calculations indicate that ^ is 
independent of rj for such a model (see, e.g., Refs. RITSi), this is strong evidence that the 
relation ^ = Dj2 holds exactly for the corresponding dynamical system. A lengthy deriva- 
tion of this relation was given in Ref. [^], based on the assumption that the pseudorandom 
disorder can be replaced by a truly random one (this gives the usual Lloyd model ||T9[|). 



The paper is organized as follows. In Sec. II, we discuss some basic aspects of the QE 
spectrum at QAR, and consider special cases of QAR occurring in ordinary KR systems. 
In Sec. Ill, we introduce the general class of modulated kicked rotors (MKRs), and deter- 
mine values of r], r^o, where QAR of period p = 1 occurs for arbitrary kicking potentials 
[more general cases of QAR, of periods p = 1 and p = 2, are considered in Appendix A; 
in Appendix B, we study QAR of arbitrary period for integrable versions of the MKRs]. 
We show, on the basis of a self-consistent approach, that for infinitesimal rj — rjQ asymptotic 
exponential localization takes place, with a pure-point QE spectrum. This spectrum and 
the QE states are determined from an effective Hamiltonian with a periodic potential, and 
the localization length ^0 is fixed by the analytical properties of this potential. If rjo = 0, 
the effective Hamiltonian turns out to be precisely an integrable limit (T — > 0) of the clas- 
sical MKR Hamiltonian. In Sec. IV, we consider cases for which exact and closed results 
concerning QAR-localization (e.g., .^0) and associated QE spectrum can be obtained. Using 
these results, we provide strong numerical evidence that QAR-localization is an excellent 
approximation to DL sufficiently close to QAR. In Sec. V, we show how MKR dynami- 
cal problems can be mapped, in general, into multi-channel tight-binding models ||20| with 



pseudorandom disorder. In Sec. VI, we study the transition from QAR-localization to DL 
in a semiclassical regime for a simple MKR system equivalent to the KR. This study is per- 
formed by considering the minimal Lyapunov exponents 'Jn for successive truncations of the 
corresponding pseudorandom tight-binding model. The exact relation = D/2 is derived 
for a nearest-neighbor pseudorandom Lloyd model. Conclusions are presented in Sec. VII. 



Some of our results have been briefly reported in Refs. 21,22 



II. QE SPECTRUM AT QAR. EXAMPLES OF QAR IN KR SYSTEMS 

An immediate consequence of is that the spectrum of U consists precisely of p eigen- 
values exp{—iui), I = 0, p — I, where the quasienergies ui are given by 

/3 + 2ttI 

uji = . (3) 

p 

Since the QE spectrum is finite, each quasienergy @ must be infinitely degenerate. An 
infinite set of QE states associated with QE level I is obtained by applying the corresponding 
projection operator for the cyclic group {e*''^/^f/*}s=o,...,p-i to all the states \1/ in the Hilbert 
space: 

iP-i 

V^z = -^e^^'^'f/"* . (4) 

P s=0 

We recall here that in the case of quantum resonance the QE spectrum consists of a fi- 
nite number of bands . The finite width of each of these bands leads to ballistic motion 
(quadratic increase of the energy expectation- value with time). In the QAR case, on the 
other hand, one has the diametrically opposite phenomenon of periodic recurrences. This 
phenomenon has nothing to do with localized QE states, since the infinite basis of states 
for QE level I can be chosen, of course, either localized or extended by properly choosing the 
state The periodic recurrences may be explained by saying that each of the p infinitely 
degenerate levels in the QAR case is the extreme limit case of a quantum-resonance band 
of zero width. This point of view will become clearer by the following examples. 
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A first case of QAR was noticed by Izrailev and Sliepelyansky in the KR. Consider 
tlie general KR Hamiltonian 

H=Yi+ kV{e)AT{t) , (5) 

where L is the angular momentum, / is the moment of inertia, is a parameter, and V{6) 
is a general periodic function of the angle 6. The evolution operator for (|^), from t = —0 to 
t = T-0, is 



where n = L/h = —id/d6, t = hT/2I, and k = k/h. Quantum resonances occur, in general, 
for rational values of = r/27r [Q. Consider, however, the special case of ?7 = 1/2. Using 
the relation 

which is easily established by comparing the Fourier expansions of V{6) and V{6 + tt), one 
finds in this case that 

t/^ = exp{-ik[V{e) + V{e + tt)]} . 

Thus the condition (D for QAR is satisfied with p = 2 if V{e) + V{e + vr) = (3/k identically, 
for some f3. This implies that V{6) must have the general Fourier expansion 



oo 



no) = ^ + E ^2s+ie^(^^+^)^ . (8) 

s=-oo 

This is, of course, the case for the standard potential V{9) = cos{9), considered in 
Ref. IQ. According to (|), the QE spectrum consists of two infinitely degenerate levels, 
uj = (3/2, (3/2 + 71. By "switching on" even-harmonic components V2s in @, the infinite 
degeneracy is removed, and the two levels broaden into two bands, corresponding to the 
generic spectrum of the 1/2 quantum resonance. 
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More general results can be obtained for the "linear" version of (|5|) 123] , which is, however, 



integrable 



r 



H=-L + kV{e)^T{t) , 
where r is now some dimensionless parameter. The evolution operator for (|^) is 

JJ — Q-iT-n^-ikV{e) 

The pth power of U in (|r^) can be easily given in closed form: 



(9) 



(10) 



UP = exp 



s=l 



exp{—ipTn) . 



(11) 



Eq. d^) is now satisfied if and only if I]s=i ~ ■5'^) = P/k and rj = t/2'k = m/p, for 
relatively prime integers m and p. The latter condition (rational rj) is precisely the general 
condition for quantum resonance P,|25[]. The former condition gives, however, the opposite 
phenomenon, i.e., the QAR. It is easy to see that this condition is satisfied only if the Fourier 
coefficients f„ of V{6) satisfy 



vo = -r, 
pk 



-'sp 



(s^O) . 



In fact, Vsp {s 7^ 0) are precisely the Fourier coefficients which contribute to the width of a 
QE band in the case of m/p quantum resonance PB . 



These examples show that in ordinary KR systems QAR may occur only if V{6) satisfies 
some restrictive conditions. In the next section we shall introduce systems in which QAR 
occurs for arbitrary kicking potentials, at some rj = rjQ. 



III. QAR IN MODULATED KR SYSTEMS AND 
ASYMPTOTIC EXPONENTIAL LOCALIZATION 

We define the general modulated kicked rotor (MKR) by the Hamiltonian 



8 



t2 M-1 

H = - + kv{e)Y.c,^T(t-t,) , (12) 

j=0 

where V{6) is an analytic function of 6*, and, for j = 0, M — 1, cj are arbitrary coefficients 
and 

0<tj < tj+i <T , to = , tM = T . 

The Hamiltonian (^) has the general form (|T|) with /(t) = C{t)Y,Y=^^ ^T{t — tj), where 
C{t) is a periodic function with period T, satisfying the M conditions C{tj) = cj. Thus 
(|l^) may be viewed as a generalized KR with M kicks at arbitrary times tj within the basic 
period, and modulated by the function C(t). The classical map for ([T^) is given by 



Ls+i = Ls- kcjV'{6s) , 

(13) 

Os+l = 9s + [{ts+l —ts)/I]Ls+i , 

where the integer s is uniquely decomposed as s = rM + j (r integer and j = 0, M — 1), 
= rMT + tj, Ls = Lit = - 0), and 6^ = 6{t = - 0). In general, the system (0) with 
(p^) is classically nonintegrable, and exhibits the transition from local to global chaos when 
k is increased, as in the ordinary KR case see an example in Fig. The simple case 
of M = 2, with Co = — Ci = 1 and ti = T/2 (the "two-sided" KR) was studied in detail in 
Refs. I^T],^. This case may already be considered as an approximation of sinusoidal driving 
potentials corresponding to ac electromagnetic fields Better approximations should be 



achieved by using the Hamiltonian (p!2D , with properly chosen coefficients Cj. The study in 



Refs. 12T|J2^ will now be extended to the general case of (p!2D 



The evolution operator for (p!2[) , from t = — to t = T — 0, is 



M-1 

U= II exp{-iTjn'^)exp[-icjkV{9)] , (14) 
i=o 



where, for j = 0, M — 1, 



= . (15) 
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and the factors under the product sign in (0) are arranged from right to left in order of 
increasing j. Now, when all the quantities (|I5|) are integer multiples of 27r, i.e., tj = 2-Kmj, 
one has exp(— zTj-ra^) = 1 identically. Then, if the coefficients Cj satisfy the condition 

M-l 

E = , (16) 

we find that f/ = 1 in (p^, corresponding to a simple case of fundamental QAR (p = 1). 
For simplicity, we shall restrict ourselves in what follows to this case, characterized by the 
conditions T.j = 271771 j and (|TB]). More general cases, for both periods p = 1 and p = 2, are 
considered in Appendix A, where it is conjectured that QAR with period p > 2 do not exist 
for the MKR Hamiltonian 



When U = 1, the QE spectrum consists just of a single, infinitely degenerate level. 
The natural question is then precisely how this infinite degeneracy is removed by slightly 
perturbing 7]j near their integer values m^. For definiteness, the perturbation of tj in (|15D 
will be made by perturbing / near I = Iq, leaving tj fixed. Denoting by e the corresponding 
perturbation in r = hT/2I, the perturbation in Tj is given by 



= '-"^e ■ (17) 



Using the operator identity [^] 



e^Be-^ = B+[A, B] + ^[A, [A, B]] + ^[A, [A, [A, B]]] + ■ ■ • , 

and formally expanding the operators exp{—iTjfP) = exp{—iejfi^) in powers of e^, we find, 
to first order in e^, 

M-l 

U^l-J2ej {m^ - d,k[27V\9)n + V"{9)] + idyV'\9)} , (18) 

3=0 

where dj = J2i=o(^s, and the prime on V denotes differentiation with respect to 9. Using 
(|T7|), the expression in ( |T8|) can be written, to first order in e, as exp(— ieGi), where 

Gi = [n - kdV'{9)] ^ + klffV'\9) (19) 
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and keff = kAd. Here d and Ad are, respectively, the average and standard deviation of dj 
with "probabihty distribution" {tj+i — tj)/T: 



d= J2 ' d 



j=0 



T 



{Ad) 



M-l ^ ^ 



E 



T 



' d] - (? 



(20) 



Using ( pUD and the definition of (ij, it is easy to show that /Cg// (or Ad) is invariant under 



cychc permutations of the sequence c,- p8 



Assuming for the moment the vahdity of the expansion above in powers of e (see the 
discussion below), the QE states in the limit of infinitesimal e are precisely the eigenstates 
of Gi, 

Gi^ = , (21) 

with quasienergies u = eg (e ^ 0). Performing on (^) the gauge transformation 

(f = exp [-ikdV{e)] tfj , (22) 

we obtain for ip, using (0), the eigenvalue equation 

-^ + k!jfr'{e)^ = gcp. (23) 

We thus see that the QE problem for infinitesimal e is just that of a Schrodinger equation 
(p3|) with a periodic potential. The spectrum g then has a band structure, but because of 
the periodic boundary condition if{2TT) = (f{0), only the level with zero quasimomentum 
is picked out from each band. This gives, in general, a point spectrum. Now, being the 
solution of the linear differential equation (pSf), ^p{0) is analytic at least in the domain of 
analyticity of V'{6) ||2^. Let 7 be the smallest distance of a singularity of V'{9) from the real 
^-axis. Then the Fourier- series expansion of ^p{0) will converge at least within an infinite 



horizontal strip of width 27, symmetrically positioned around the real 6'-axis It follows 
that the Fourier coefficients of ip and ip in (^) decay asymptotically at least as exp(— 7|n|). 
This means that in the immediate vicinity of QAR there takes place asymptotic exponential 
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localization in the angular momentum nh, with locahzation length not larger than I/7. 
In general, .^0 is determined entirely by the analytical properties of V'{6) (see examples in 
next section). 



This exponential "QAR- localization" in L, following from Eq. (P5|), justifies a posteriori 
the expansion above in powers of e. In fact, the general expansion for U in ([l^ ) can be 
formally written as exp(— iG), G = Yl,JLi^''Gj. Here the Hermitian operators Gj are poly- 
nomials in n and derivatives of V{9) of order not larger than 2j [the leading operator, Gi, 
is given by (0)]. Thus, the highest power of h contributed by Gj appears in this expansion 
as {efpy . This means that by choosing e < n~^^, where Umax ^ I/7, the eigenstates of 
U (i.e., the QE states) should be very close to those of Gi, at least within the localization 
domain. In the limit of infinitesimal e, the QE states should coincide with the eigenstates 
of Gi. The derivation of Eq. (p3|) appears then to be self-consistent. 



We now show that the effective Hamiltonian (|T^) has a classical counterpart in the 
limit of very small values of the quantity T/I = 2T/h. This limit corresponds to the 
case of infinitesimal values of tj in (|r^), i.e., infinitesimally close to the special QAR point 
rj = t/27t = {nij = for all j). Consider the Mth iteration of the classical map (|13|), 
giving the map {Ls, Og) {Ls+m, (^s+m)- Taking carefully the limit T ^ in this map 
at fixed (tj+i — tj)/T and using the condition (|16]), we obtain, after a straightforward but 
tedious calculation, the Hamilton equations 



dL _ dHeff de _ dHeff 
It ~ W ' It ~ dL 



(24) 



where 



Heff = Yi^^ = Yj\L- kdV'iO)] ' + ^^^V'\e) . (25) 
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Eqs. ( P^ and ( P^D show that the general MKR Hamiltonian (|T2|), with coefficients Cj sat- 
isfying (|16D, is integrable in the limit T — 0, as it reduces precisely to the ID effective 
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Hamiltonian (pS]). The latter is essentially the QAR effective Hamiltonian (|T9]), and, af- 
ter the canonical transformation L' = L — kdV'{d) [analogous to the gauge transformation 
(p^)], it becomes essentially the Schrodinger Hamiltonian in p3|). Thus, QAR-localization 
in the infinitesimal vicinity of r/ = is associated with a classically integrable limit. As T 
is increased from 0, keeping the quantities (tj+i — tj)/T fixed at some rational values irtj/m 
(m = Z^jlo^ QAR-localization for = will repeat periodically in the infinitesimal 

vicinity oi r] = rm, for all integers r. For these values of i], which are equivalent to rj = 
but correspond to nonintegrable systems in a strong quantum regime, the QAR-localization 
is only a "reflection" of the classically integrable limit T — 0. In Appendix A, we show that 
QARs of periods p = 1 and p = 2 can occur, in general, if tj in (|15D is an odd multiple of 
71. Such values of tj are not equivalent to tj = 0, since tj mod 27r ^ 0. In this case, the 
QAR-localization is not even a reflection of a classically integrable limit. 

It is important to notice that the limit T — > (or r/ ^ 0) at fixed k (or k) is not a 
semiclassical limit. In fact, if the quantities (tj+i — tj)/T are kept fixed and the coordinate 
transformation L' = {T/I)L is performed in the map ([TB|), it becomes clear that the classical 
dynamics depends only on the parameter K = {T/I)k = 2Tk. The semiclassical limit is 
then T] at fixed K, not at fixed k. However, at fixed A; ^ 1, small values of r] such that 
K may be viewed as corresponding to a semiclassical regime of almost integrability. 

IV. EXACTLY-SOLVABLE CASES 

In this section we consider cases of potentials V{6) for which the QE problem in the 
infinitesimal vicinity of QAR [Eq. (|23|)] can be solved in closed form, or at least an explicit 
expression can be obtained for the QAR-localization length ,^o- Using these exact results, 
we shall provide strong numerical evidence that the QE spectrum and localization features 
sufficiently close to QAR are well accounted for by the QAR effective Hamiltonian (|T^). As 
shown in the next section, the MKR dynamical problem can be mapped into pseudorandom 
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tight-binding models, so that dynamical localization (DL) is expected to occur if 77 = r/27r 
is sufficiently irrational. If, in addition, rj — tjq is small enough, this DL should look similar 
to QAR-localization. 



Our ffist example is the standard potential V{9) = cos (6*), for which Eq. (^) reduces to 
the Mathieu equation ||29|j30( 



y" + [a-2q cos(2^)] y = , (26) 

where y = (p, a = g — {kAd)'^/2, and q = —{kAdy/4. The problem is then exactly solved 
in terms of the periodic Mathieu functions y = cCriO, q) (symmetric) and y = ser{0, q) 
(anti-symmetric), with corresponding eigenvalues a = ar{q) and a = br{q)- Explicit expres- 
sions for these functions and eigenvalues, as well as a detailed discussion of their properties, 
can be found in Refs. [p9| , p0| . From Eq. ( ^2]) the Fourier coefficients ipn and ?/„ of ip and 



y = Lp, respectively, are related by ipn = Hji^ Jj{.k(^)yn-ji where Jj{kd) is a Bessel function. 
Since the dominant decay rate of both J„ and ?/„ with n is like n"" [^, this is also the 



dominant decay rate of ipn- This strong localization in L space, faster than exponential, 
could be expected from the fact that V{6) = cos{6) is an entire function (analyticity-strip 
width 27 = 00), so that the asymptotic localization length ^0 = I/7 = 0. 



Let us now check to what extent the QAR effective Hamiltonian (|T^ for V{6) = cos{6) 
reproduces accurately the quantum dynamics and QE spectrum for rj — rjQ = e/27r sufficiently 
small and irrational. We have studied numerically the case of M = 3, with cq = ci = 1, 
C2 = —2, and tj = r/3 for all j (chaotic orbits for this system are shown in Fig. |l|). The 
quantum dynamics of a wave-packet initially equal to |n = 0) was investigated using a basis 
of up to 512 angular-momentum states around n = 0. The wave-packet was propagated 
in time using well-known algorithms In Fig. |^, we plot the kinetic-energy expectation 
value, Eq = (L^/2/), as a function of the "real time" t = es {s is the number of applica- 
tions of U), for k = 2.0 (fee// = kAd = -y/s/S) and several irrational values of e/2ii. We 
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observe that almost all the data falls close to the same curve, even for values of e as large 
as e = 0.28. For smaller values of fee//, the data falls much more accurately on the same 



curve (see an example for M = 2 in Ref. ||22|). This is evidence that, even for e not very 



small, the quantum dynamics over a significant time interval can be well described by the 
approximate evolution operator U'^ = exp(— iesGi) = exp(— -itGi), generated by the QAR 
effective Hamiltonian Gi. In Fig. we plot the Fourier transform Eq^u) of EQ{t) for e = 0.1 
and several values of k. The positions of the various peaks in Eo{i>) must correspond to 
values of eu equal to the spacings between QE levels. A comparison with the level spacings 
corresponding to eigenvalues of the Mathieu equation (|26|) shows excellent agreement. This 
is strong evidence that sufficiently close to QAR the QE spectrum is very well accounted 
for by the QAR operator d [i.e., Eq. (H)]. 

Dynamical localization (DL) in the vicinity of QAR will now be compared with 
QAR- localization (we assume in what follows that rjQ = 0). For this purpose, it will be 
sufficient to consider the localization of steady-state probability distributions {\'ip{n)\'^) over 
angular momentum nh. Such a distribution is calculated by propagating an initial angular- 
momentum state \n = 0) for a sufficiently long time interval, and averaging then the results 
for |'?/'(n)p over a subinterval at the end of this interval. All our numerical calculations in 
what follows have been restricted to the simple case of the M = 2 MKR with cq = —C\ = 1 
and ti = T/2 (the "two-sided" KR |pT| , |22| ). Fig. 4 shows a semi-logarithmic plot of {\ilj{n)\'^) 
for V{6) = cos{6), k = 10, and three irrational values of 77 = r/27r, r = 10"^, 10^'^, 10^^. 
For the first two values of r, corresponding to close vicinity of QAR, {\ilj{n)\'^) appears to 
decay faster than exponential, in accord with the solutions of the Mathieu equation (p6|) (i.e., 
QAR- localization). On the other hand, the localization for r = 10^^ looks quite different 
from QAR- localization; it appears, in fact, to be almost exponential. These results can be 
explained as follows. The classical nonintegrability parameter is given by = (T/I)k = 2rA; 
(see previous section). Thus, the first two values of r correspond to a semiclassical regime 



(small r) of almost integrability (small K), and the observed DL [|T^] is then similar to 
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QAR-localization, as predicted from the integrable effective Hamiltonian ([TOD (however, a 
basic difference between these two locahzations is expected, as pointed out in Sec. VI). 
For r = 10^^, on the other hand, K = 2, corresponding to a semiclassical regime (small 
r) of global chaotic diffusion in the M = 2 MKR As in the KR case the lo- 

calization observed here should be characterized by a DL length ^ completely determined 
by the chaotic-diffusion coefficient D, ^ ^ D/2 (see Sec. VI). This is in contrast with 
QAR-localization, which is determined by the analytical properties of the potential. The 
transition from QAR-localization to DL, as r is varied in a semiclassical regime at fixed k, 
will be studied in more detail in Sec. VI. 



Our second example is the potential 

V{e) = Aarctan[Kcos(^) - kq] , (27) 

where A, k, and kq are some constants. The QAR-localization length for (0) can be 
determined exactly as follows. The function V'{6) assumes simple poles 6q satisfying the 
equation kcos(^o) ~ ^^o = i^- The distance of any of these poles from the real 6'-axis is 
7 = |Im(^o)|) and we easily find that 

2k cosh(7) = [1 + {ko + «:) + [1 + («o - «;)T^' • (28) 

Consider now the Fourier-series expansion for the solutions <f{0) of Eq. (^). We claim that 
in the case of ( pT]) the Fourier coefficients <y9„ must decay asymptotically as (fn ^ 6xp(— 7|n|). 
The QAR-localization length is then = 1/7- To show this, we observe that the simple 
poles 6q correspond to regular singularities of order 2 of Eq. (pSf). The exponents pi 
and p2 for these singularities are easily determined from the quadratic "indicial" equation 
p9| for (pBI): pi = p2 = 1/2. Since the exponents are equal, the general solution of ( |23D 
around 9 = 9q assumes the form ||2^ 



m = {9- 9oy/^Ri{9 - 9o) + R2{9 - 9o)[W + h \og{9 - 9^)]} , (29) 
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where Ri{9) and R2{0) are analytic (can be expressed as Taylor expansions) around 9 = Oq, 
and bi and 62 are arbitrary constants. It follows from (|29| ) that all the derivatives of <f{0) 
diverge at = ^^o- Let us now continue the Fourier- series expansion for ip{6) into the com- 
plex 6 plane. Defining the complex variable z = e*^, one gets a Laurent expansion in z 
which converges at least in a "ring" excluding the singularities of Eq. ( P^ pP| , i.e., for 
|Im(^)| < Im(^^o) = 7- However, since all the derivatives of (^91) diverge as — 6q, this must 
be also the case for the derivatives of the Laurent expansion. By a simple consideration of 
the latter derivatives, the desired relation ipn x exp(— 7|n|) is obtained. 

Figs. 5 and 6 show semi-logarithmic plots of (IV'(^)P) for the M = 2 MKR (defined as 
above) with the potential (pTj) {A = 1, k, = 1, and kq = 0) and for several values of k and 
r. The straight line (dotted) in both figures has slope —27, where 7 is determined from Eq. 
(p8|). For T sufficiently small, this slope is expected to be close to the asymptotic rate of 
exponential decay of (|^/'(n)p). We see that this is indeed the case whenever the classical 
parameter K = 2rA; is small enough, K < 1, corresponding to almost integrability or local 
chaos. In particular, for A; = 2 (Fig. 5) the decay rate appears to be equal to 27 for all 
three values of r. In fact, for k = 2 and A = 1 the potential ( P?] ) leads to a nearest-neighbor 
pseudorandom tight-binding model whose DL length ^ appears to be independent of r (see 
Sec. VI and Fig. 8). On the other hand, for k = 10 (Fig. 6) the DL length is quite sensitive 
to the value of r. In the case of r = 10^^, corresponding to a semiclassical regime of global 
chaotic diffusion {K = 2), the observed DL looks quite different from QAR- localization, 
with an asymptotic DL length ^ ^ D/2 ^ C,o- This case is similar to that of r = 10~^ in 
Fig. 4. 



V. MULTI-CHANNEL PSEUDORANDOM 
TIGHT-BINDING MODELS 

We now show how the MKR dynamical problem with ( |T2|) can be mapped into a pseu- 
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dorandom tight-binding model, in analogy to the ordinary KR case ||^. For simplicity, we 
shall assume that the quantities (|T5|) are all equal, tj = t/M for all j (i.e., the kicks are 
equidistant in time). Let uf{6), j = 0, M — 1, denote a QE state with quasienergy cu at 
time t = tj ± 0. The following relations hold: 

u+{9) =exp[-iCjkV{9)]u- , (30) 



u 



e— ^/^^«t (0 < J < M - 1) , = e^'^—^^^t.-i.^ , (31) 



where is the L-representation of uf{9). We define, in some analogy with Ref. 



(32) 



so that Wj{e) = -tan[cjkV{e)/2]. Simple manipulations of Rels. (|30|), (|TD, and (H) yield 
a system of M equations for the L-representation Uj^n of Mj(6'), j = 0, M — 1: 

r \ r / 

where UM,n = UQ,ni and Wj^n is the L- representation of Wj{6). Unless otherwise specified, 
the index r in (^) runs over all the integers. We now introduce the Fourier transforms Ug^n 
and W s,n of Mj> and Wj^n, respectively, in the variable j: 

M-l 



1 

M 



^.,n = - E u,,ne^^'''"' , (34) 
i=o 



1 M-l 

Ws,n = V7 E W^.,ne^-^-^/^ . (35) 
3=0 



Using the expressions (plf) and (pSj) in (|33|), we obtain, after simple algebraic manipulations. 



Tlf^Us,n + E W^O,n-rMs,r + E E W^s-s',n-rMs',r = -Eu^.n , (36) 



where, for s = 0, M — 1, 



T^") = - tan[(™2 - 2713 - u)/2M] , (37) 



the index s' takes all the integer values s' = 0, M — 1 with the exception of s' = s, and 



E = —Wofl- Eqs. (|36|) describe a tight-binding model of an M-channel strip [|20|. The 
on-site potential in channel s is given by Tlf\ while the hopping constants within a channel 
are W^o,n- The hopping constants from channel s to channel s' ^ s are given by W s-s',n- 

A particularly interesting case arises when Cj assumes only the values 0, ±c for some 
constant c, which will be chosen, without loss of generality, equal to 1. From Rel. (|35|) and 
the definition Wj{9) = -ta.n[cjkV{9)/2], it follows then that 

1 M-l 

Ws,n = C(s)iy„ = V7 E Cje'^'^'^^'Wn , (38) 

j=o 

where c(s) is defined by ( |38D and Wn is the L-representation of W{9) = — tan[A;V"(^^)/2]. Rels. 
(^) and (|l^) imply now that the hopping constants M^cn within a channel are identically 
zero, including, of course, E = — VTo.o- In this singular case, the model (|36D looses much of its 
physical meaning. In Appendix C, we consider in some detail nearly-singular cases of (p6[), 
for which the hopping constants Wo,n are nonzero but small. Notice that the M = 2 case 
is always singular, since cq = — ci from (|16|). In this case, however, a physically meaningful 
two-channel model can be derived directly from Eqs. (B^), without the need of the 



Fourier transforms (|3^) and (p5|). Simple manipulations of Eqs. (|33D give in this case: 

TnUo^n + SnUi^n + Z^r^O ^n-r^Cr = EUq^u i ^^^^ 

where T„ = cot(Cri,), Sn = —1/ sin(^„), (n = — ^)/2, and E = —Wq. The on-site 
potential and hopping constants within each channel are, respectively, Tn and Wn, while Sn 
are the coupling constants between the channels. As shown in Sec. VI, the M = 2 MKR 
is essentially equivalent to an ordinary KR if the potential satisfies V{6 + tt) = —V{6). 
In this case, the dynamical problem can be conveniently approached using the well known 
tight-binding models for the KR 
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For irrational r] = r/27r, T^"^) in (|37|) [or T„ and S'„ in (|39|) ] is a pseudorandom sequence 
PJTEIl which, by arguments similar to those used in the KR case , may lead to Anderson- 
like localization of the eigenstates of (|3^) [or of (p9|)]. This is the dynamical localization (DL) 



T7| of the corresponding QE states in angular momentum. Another source of randomness, 



which does not involve r] but may contribute significantly to DL, is the dependence of the 
hopping constants Ws,n on the distance s between channels, especially in the limit of very 
large M. This may reflect a possible randomness of the sequence Cj in the "time index" j. 



Clearly, the pseudorandomness is not defined in the infinitesimal vicinity of QAR (in- 
finitesimal e = r mod 27r), so that QAR-localization is not DL, strictly speaking. However, 
when approaching a QAR point t]o, the pseudorandomness is guaranteed by choosing e/27r 
in a sequence of "strong" irrationals ei/2'K, 1 = 1, 2, converging to [e.g., ei = 27r/ (/ + g), 
where g is the golden mean]. It is then quite possible that for e = ei the model (|36D [or (p9|)] 
will have exponentially- localized eigenstates. In fact, the numerical evidence presented in 
the previous section (see also next section) strongly indicates that this is indeed the case for 
/ sufficiently large. Then, provided the QAR condition (|16D is satisfied, the quasienergy uj 
in (^^ should be given approximately hj u = eg, where g is an eigenvalue of the leading 
operator Gi in (p^. Moreover, in the M = 2 case of (^), one can easily determine, using 
(pop, (p2D, and (^2]), the accurate relation expected between the solutions of Eqs. ( |23D and 
(§|) for small e: Uo{e) ^ cos[kV{e)/2]^{e). 



For rational values of rj = m/p {m and p are relatively prime integers), the on-site po- 
tential Tjf^ in (|37D [or T„ in (0)] is periodic in n with period p. Then, by straightforward 
application of Bloch theorem, one can easily show that the QE spectrum uj is, in general, 
absolutely continuous, consisting of p bands. This corresponds to usual quantum resonance 
1^. However, when the QAR conditions are satisfied, either for p = 1 or p = 2 [see (|1^) and 
Appendix A], a QE band "collapses" to an infinitely degenerate level. 
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VI. TRANSITION FROM QAR-LOCALIZATION TO 
DYNAMICAL LOCALIZATION IN A SEMICLASSICAL REGIME 

In this section, we study the dependence of the DL length ^ on the distance r] — r]o from 
the QAR point rjo = 0. This study will be performed in a natural way at fixed k, so that 
r^o = is associated with a classically integrable limit (see Sec. III). If k is sufficiently 
large, the classical parameter K = 2rA; = Anrjk may correspond to a case of global chaotic 
diffusion already for small values of 77 < 1, i.e., in the semiclassical regime. The approximate 
relation ^ ^ D /2 where D is the diffusion coefficient [see definition (^9[) below], should 
then hold. The transition from QAR-localization to DL in the semiclassical regime, as i] is 
increased from 0, can be understood in a most illuminating way if the problem is approached 
by using a suitable pseudorandom tight-binding model, equivalent to the dynamical system. 
This approach is also most convenient for an accurate calculation of ^. 

For simplicity, we shall restrict ourselves to the M = 2 MKR (cq = — ci = 1, ti = T/2) 
with potentials satisfying ¥{6 + 71) = —V{6) (as in the numerical examples in Figs. 4-6). We 
denote by U2kr{t) the evolution operator for this system at a given value of r. Similarly, 
we denote by Ukr{j) the evolution operator (^) for the ordinary KR. Using Rels. (|T^) and 
(0) with V{9 + -k) = —V{6), it is easy to derive the following exact relation 

U],j,{7r + e/2) = U2KRie) . (40) 

Rel. (^) means that the quantum dynamics of the M = 2 MKR at distance e from the 
QAR point r = is essentially equivalent to that of the ordinary KR at distance e/2 from 
r = TT. The latter value of r is precisely a period-two QAR point for the KR (see Sec. II and 
Appendix A). This equivalence between the two systems enables one to study the M = 2 
MKR using the simple tight-binding models for the KR instead of the two-channel 

model (|5|). 
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Our basic potential is the standard one, V{6) = cos{6). A convenient tight-binding 
model for the KR with this potential was proposed by Shepelyansky [0: 

oo 

J2 Mk/2)sm[{Tn^ -uj + nr)/2]un+r = , (41) 

r=— oo 

where Jr{k/2) is the Bessel function, u is the quasienergy, and is related to the angular- 
momentum representation of the QE states. Using Rel. (^) in (^T]), we obtain the corre- 
sponding tight-binding model for the M = 2 MKR: 

oo 

Mk / 2) sm{[{27r + T)n^ -tu + 27ir]/A}un+r = , (42) 

r=— oo 

where all the quantities r, u, and Un now refer to the M = 2 MKR. Let us briefly recall 
how the asymptotic DL length C, can be calculated from such tight-binding models 
The hopping constants in (0), given by the Bessel function Jr{k/2), decay faster than 
exponentially for |r| > k/2: 

{A/nky/^ cos[A;/2 -Trr/2- n/A) for |r| < k/2, 
Jr{k/2) ~ { (43) 
'l/27r|r|)i/2(g^/4^)|r| |^| ^ ^/2. 



It is therefore reasonable to approximate (|42|) by restricting r to the finite range |r| < A^, 
for sufficiently large A^. This truncated form of Eq. (^21) can be easily written as a transfer- 
matrix problem 

as+i = Tsas , (44) 

where as is the 2A^-dimensional vector with components cr^''^ = Ug-r, r = —N + 1,...,A^, 
and Tg is a 2N x 2N symplectic matrix. One may interpret (^) as a map describing a 



Hamiltonian dynamical system with degrees of freedom |T3. The vector ctq is mapped 
into 0"„, for arbitrary n > 0, by the product matrix 

An = r„_ir„_2 ■ ■ ■ To . (45) 

Since the matrix ( ^5]) is, obviously, symplectic, its eigenvalues X{n) come in N reciprocal 
pairs [Ar-(n), A7"'^(n)], r = 1,...,A^, and we can always assume the ordering 1 < |Ai(?2)| < 
1^2(^)1 < • • • |Ajv(r;,)|. The minimal Lyapunov exponent for the map (0), 
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77V = lim -ln|AiH| , (46) 

n^oo n 

determines then an A^th-order approximation ^tv = l/7Af to the asymptotic DL length ^. 
Since an accurate calculation of becomes extremely time-consuming as is increased, 
important questions are how fast C,n converges to its limit value ^ as is increased, and 
whether has a well-defined quantum-dynamical meaning per se. 



We have therefore studied as a function of both N and t] near QAR. The value of k 
was fixed at k = 10, and several values of r were considered in the interval 10^^ < t < 1. 



We have calculated 77V for A^ < 15, using the well known method |]7|,|8|jl8[| for determining 
the Lyapunov spectra of products of matrices such as (^Sf). The method is based on direct 
application of the map (|^) a large number n = rimax of times, such that for n > rimax the 
matrices r„ can be considered as random. This randomness should be realized to some ex- 
tent by the pseudorandom term in (^) if r[(?2+l)^ — n^]/4 ~ 27r, or n = Umax ~ ^ti/t. 
We have used rimax = 10"^ for r = 10^^, 10^^, and rimax = 10^ for the larger values of r. 
We have checked that this choice of rimax yields well-converged results by calculating and 
comparing 7Ar for different values of n < rimax, e.g., n = 10^, 2 ■ 10^, 10^ for Umax = 10^. 
The final results are shown in Fig. 7 using a semi- logarithmic plot. We observe that 7Ar 
clearly decreases in the interval A^ < 3 for all values of r. For larger values of A^ and for very 
small T (immediate vicinity of QAR), 77V appears to increase without bounds with A^. For 
r > 0.1, 7Ar appears to "saturate" around some limit value which decreases as r is increased. 
For r ~ 1, the saturation occurs immediately after A^ > 3. 



The apparently unbounded increase of 'Jn with A^ for very small r is consistent with the 
fact that the asymptotic QAR-localization length for V{6) = cos{6) is ^0 = (analyticity- 
strip width 27 = 2700 = 00). We now show that 77V has a well-defined quantum-dynamical 
meaning also for finite A^. This is because the truncated version of the tight-binding model 
(^) turns out to be exactly equivalent to the dynamical problem (KR or M = 2 MKR) 
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with a potential Vn{0) replacing V{9) = cos(^). Then = '^Iin is the DL length for this 
potential, reducing to the QAR-localization length ^^vo in the limit r ^ 0. For very small 
r (e.g., r = 10~^ and r = 10~^ in Fig. 7), should be an excellent approximation to ^at^o- 
The potential Vn{0) can be easily determined from the (7-function approach of Shepelyansky 

N 

WN{e) = gN{e)eM-^kVN{e)/2] = Jr{k/2)d<'-^/^'^ , (47) 

r=-N 

where gNifi) is some real function. Solving Eq. (^7|) for Vn{0), we obtain 



2 

^Af(^) = — — arctan 
k 



E.=_w^r(fc/2) sin(r0 - r7r/2) 



(48) 



Xr=-N Mk/2) cos{re - r7r/2)_ 
Using J_r{k/2) = {—iy'Jr{k/2), it is easy to see that (^) satisfies Vj\f{9 + vr) = —V]\f{9), 
so that Rel. (|0|) between the KR and the M = 2 MKR holds for the potential ¥^{9). For 



= 1, (|48|) reduces to (pTj) with = (see below). The QAR-localization length for (^8]) 
can be determined as in the case of the potential (pTj). It is easily shown that Vl^{9) assumes 
simple poles 9o satisfying the equation 

N 

WNi9o) = E Jrik/2)e'^^'^'~-/'^ = , 

r=-N 

or gN{9o) = 0, since the exponential function in (^) can never vanish. The QAR-localization 
length for Vn{9) is given by ^jvq = I/Tm'^ where 7^"* is the smallest distance of a pole 9^ 
from the real 6'-axis. It is also the half- width of the strip of analyticity of Vn{9) (Fourier- 
series representation). Since 7 = 00 for V{9) = cos{9), MniN^oolN^ = ^5 so that 7^^ must 
generally increase with A^, as is clear from Fig. 7 for very small r. 



We thus see that in the limit r — > 7Ar is both the minimal Lyapunov exponent of a 
2A^ X 2A^ symplectic matrix and the half- width of analyticity of the potential Vn{9) for the 
truncated problem. As A^ increases, the region of analyticity of Vn{9) increases without 
bounds, and the QAR-localization length ^7v,o ~^ 0. With this in mind, we now consider the 
cases where ■jn appears to saturate for A^ not too large in Fig. 7. These cases correspond 
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to regimes of local {K < 1) or global {K > 1) chaos. As in the KR case [|^,|^, we expect 
that in the global-chaos regime the saturated value of the DL length ^f^f = l/'j^ is given 
approximately by ^ ~ D/2, where D is the classical chaotic diffusion coefficient. In complete 
analogy with Refs. 0J^, we shall define D for the MKR map (|13|) as follows: 

D = hm ^^^rM- L,r) ^49^ 

where ( ) denotes average over an ensemble of initial conditions {{Lq, 6o)}- In a strong-chaos 
regime, we may assume, as usual, that the angles dg are independent and random, giving 
vanishing force-force correlations, {V'{6s)V'{9s')) = 0. Using then (|T^) in (49), we obtain 
the following expression for D: 



k\ 1 

D = — lim - 

r— ►oo y 



r-1 M-1 
r'=0 j=0 



e{c') f\'\9)^, (50) 



where (c^) = J^fLo^cj/M, and, as in Ref. |^, we use units such that ^ = 1, giving k = k. 
For the M = 2 MKR (c^) = 1, and, with A; = 10 and ¥{6) = cos{e), we find from (0) that 
D = 50. The minimal saturated value of 77V in Fig. 7 is ~ 0.035, so that the maximal 
DL length is ^ ^ 28 ~ D/2. We have checked that the approximated relation ^ D/2 holds 
also for other values of k. 



We observe that the minimal saturated value of ~ 2/D is attained immediately after 
the decrease of 77V in the interval < 3. This behavior was observed also for other values 
of k, with in the interval A^ < k/2. The decrease in 77V may be due to the fact that for 
\r\ < k/2 the hopping constants (|^) are of the same order of magnitude, so that the effect 
of each additional hopping constant is to increase the DL length. For |r| > k/2, on the 
other hand, the hopping constants ( ^3] ) decay faster than exponentially, and their effect on 
7Ar turns out to be negligible in the global-chaos regime, as shown in Fig. 7. This fact was 
used in Ref. |^ to justify the truncation of the tight-binding model for the KR at N ^ k/2. 

We see in Fig. 7, however, that the effect of the exponentially small hopping constants 
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for \r\ > k/2 is not negligible for small values of r < 1, corresponding to local-chaos regimes 
or the vicinity of QAR. In sharp contrast with the case of < k/2, these small hopping 
constants have the somehow paradoxical effect to increase 'Jn, i-e. to decrease the localiza- 
tion length C,- This effect appears to continue up to a maximal value A^* of A^. For A^ > A^*, 
7Ar saturates to a value ~ jiy*. Both A^* and 77V* increase as r or /c decrease. 

While we are unable to provide at this point a quantitative explanation of these phenom- 
ena, it will become apparent from the following discussion that they are quite natural. We 
already know that the truncated tight-binding model is exactly equivalent to a dynamical 
problem with a potential Vn{9) [see p8|)]. The region of analyticity TZjy of Viy{9) increases 
with A^ for A^ sufficiently large (A^ > k/2). The strip of analyticity 71qe,n of the QE states 
has a width 27jv, where 77V depends on r and determines the rate of exponential decay of the 
QE states in angular-momentum space. In the infinitesimal vicinity of QAR, 1Zqe,n = T^n, 
so that an increase in the analyticity of the potential results in a corresponding increase 
of the analyticity of the QE states. For finite r, however, the increase of 7^ at leads to an 
increase of 'JZqe,n only up to A^ ~ A^*, where saturates to a value ■jjy*. As the strength of 
chaos is increased by increasing r (or K), 77V* decreases (the DL length ^ increases). At the 
same time, A^* decreases, so that the infiuence of the analyticity of Vn{0) on the analyticity 
of the QE states is gradually reduced. The extreme case corresponds to the global-chaos 
regime. Here the increase of TZ^ does not lead to any increase in 71qe,n, and 'Jn* is totally 
determined by the diffusion coefficient D. This case of DL in a semi classical regime of global 
chaos is thus completely different in nature from that of QAR-localization. We may now 
say that the transition between these two kinds of quantum localization takes place via a 
gradual reduction of the influence of the analyticity of the system on that of the QE states, 
as the level of chaos is increased. 

This characterization of the transition is quite natural as it has a well-known classical 
analogue. When a nonintegrability parameter is increased, the analyticity of functions rep- 
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resenting classical structures (e.g., KAM tori, periodic orbits) generally decreases. A famous 
example is the golden-mean torus in the standard map. When the parameter K approaches 
from below the critical value ~ 0.9716 for the disappearance of this torus, the width Agm 
of the strip of analyticity of functions representing the torus shrinks to zero a,s Kc — K (for K 
sufficiently close to |^. Ki K = Kc, the torus is a continuous nondifferentiable curve, 
and for K > Kc it becomes a cantorus allowing for global chaotic diffusion. The decrease 
of Agm as K is increased is analogous to the decrease of 7Ar*. Unlike Aqai, however, the 
minimal value of 7^^*, in the global-chaos regime, is not zero but only inversely proportional 
to D > 1. 



According to this picture, a saturation of '-/n around an asymptotic value ~ 'Jn* should 
occur for arbitrarily low level of chaos (arbitrarily small r or /c), but A^* and 'Jn* may be 
very large. We thus expect that, unlike QAR-localization, the DL decay can never be faster 
than exponential but it should always feature an asymptotic exponential tail. For very small 
r, this asymptotic behavior may start only at very small value of the wave function, and the 
decay may look faster than exponential above this value. This is probably what one sees in 
Fig. 4 for T = 10~^ and r = 10^'^. Here an exponential tail may be observed only much 
below the level of numerical noise (~ 10^'^'^). Because of computational limitations, we were 
not able to verify the existence of the saturation effect and the associated exponential tail 
in the DL decay for r < 0.05. 



In Fig. 8, we plot 'Jn for A^ = 12 and A^ = 1 as a function of r near QAR. While 712 
changes significantly from its maximal value (near QAR) to its minimal DL value ^ 2/D, 
7i appears to be independent of r. In fact, the small fiuctuations in 71 are of the same 
order of magnitude as the numerical error in this quantity. The case of A^ = 1 corresponds 
to a nearest-neighbor pseudorandom "Lloyd model" PJ^,P,p!5|] . The independence of the 
localization length ^ = I/7 on r for such a model has been verified numerically for 
many values of a coupling constant (e.g., k). The potential ^1(6*) in (^Sj) is a special case of 
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mi 

2 

Vi{9) = — -arctan[Kcos(0)] , (51) 

rh 

where k = —2Ji{k/2)/Jo{k/2). For A = —2/k and arbitrary hq, the potential ( pTj) leads 
to nearest-neighbor hopping constants also in the original M = 2 tight-binding model (p9[): 
Wn = except of W±i = k/2. This is a two-channel pseudorandom "Lloyd model" (see 
some generalizations of this model in Appendix C). 

In Ref. the relation ^ = D /2 was derived for the KR with the potential (|5TD by assum- 
ing that the pseudorandom disorder can be replaced by a truly random one, giving the usual 
Lloyd model [|l^. The validity of this assumption was verified in Ref. |jl5|. We now give an 
alternative and simpler "proof of this relation, based on the exaci result C,o = D/2 which is 
derived below. The relation ^ = D/2 is an immediate consequence of the latter result and 
the numerically observed independence of 71 on r. The derivation below of .^0 = -D/2 makes 
use, of course, of the fact that the KR with the potential ( |5TD exhibits QAR (of period two). 
As shown in Sec. II, this QAR can occur only if V{9 + 7t) + V{6) = constant. This condition 
is satisfied by the potential (^) only if /to = with constant = 0, which is precisely the case 
of (^l|). As shown below, the relation ^ = D/2 is valid also for the equivalent system of the 
M = 2MKR with (|3). 



To derive the relation ^0 = D/2, we first perform explicitly the integral in (|50D for the 
potential (|5TD. The final result can be easily expressed in terms of 7 in (pSl): 



D = 2{c^) 



- + ^ + o(y) 

7 2 



(52) 



where we have assumed the strong-chaos regime k ^ 1 for (pi]), corresponding to very small 
7 in (|28D. Now, the QAR-localization length ,^0 for ( F^D is ^0 = 1/7 (see Sec. IV), and 
(c^) = 1 for both the KR and the M = 2 MKR. The relation ^0 = D/2 for these systems 
follows then from Eq. (B2p in the strong-chaos regime. 
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From Fig. 7 it appears that 72 is also independent of r for the particular value of k 
chosen. The case of = 2 corresponds to a next nearest-neighbor pseudorandom model. 
Such a model was investigated by Brenner and Fishman [|I5| who found that its localization 
properties are quite different from those of its truly random counterpart. In particular, 
the localization length strongly depends on the irrational number (analogous to rj = t/2'k) 
defining the pseudorandom disorder. Fig. 7 suggests, however, that it may be possible to 
define pseudorandom tight-binding models with > 1 neighbors, whose localization length 
is almost independent of r. Such models should exhibit localization properties quite similar 
to those of their truly random counterparts. The existence of these models and related 
problems will be investigated in a future work p2| . 
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VII. CONCLUSIONS 



In this paper, the problem of quantum locahzation in 1.5 degrees of freedom ("minimal 
chaos") has been studied in a large class of systems exhibiting classical chaotic diffusion, 
the modulated kicked rotors (MKRs). These systems feature two basically different kinds 
of asymptotic exponential localization in angular-momentum space: (a) Dynamical localiza- 
tion (DL), which, as in the case of the KR and other systems P,f7|,^ |lOi|TI| , is the localization 
exhibited by pseudorandom tight-binding models for irrational values of 1] (a scaled h) ||17|| . 
(b) QAR-localization, which is the localization occurring for t] in the infinitesimal vicinity of 
points Tjo where the quantum dynamics is exactly periodic [quantum antiresonance (QAR)]. 



The existence of QAR-localization has been rigorously established in the framework of 
a self-consistent approach for both period-one and period-two QAR (see Sec. Ill and Ap- 
pendix A; it seems that QARs of period larger than two do not exist for nonintegrable 
MKRs). This approach leads to the basic equation ( p3D [or ( |A11| ) in Appendix A] for the 
QE problem in the infinitesimal vicinity of QAR. It follows from this equation that the QE 
spectrum is pure point for infinitesimal f] — f]Q, and that the asymptotic QAR-localization 
length ^0 is completely determined from the analytical properties of the potential appearing 
in the equation: ^ 1/7; where 7 is the smallest distance of a singularity of the potential 
from the real ^-axis. In many interesting cases, such as the potentials (|27D and (0), ■Co = I/7 
exactly. 



Being associated with the ID time- independent Schrodinger equation (p3|), QAR- 
localization is of an "integrable" nature. In fact, if r^o = is a QAR point, the Hamiltonian 
in (|23| ) is precisely the integrable limit (T ^ at fixed k) of the classical MKR Hamil- 
tonian (p^). On the other hand, DL is associated with the difference equations ( |36D [or 
(p9|)] for the pseudorandom tight-binding models. These equations depend not only on the 
dimensionless kicking parameter k [essentially the only parameter in (^)], but also on 77, 
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which determines the pseudorandom disorder. This pseudorandomness, which is absent in 
the angular-momentum representation of (|23|), introduces "nonintegrable" features in DL. 



Thus, while the asymptotic QAR-localization length ,^0 is completely determined by the 
analytical properties of the potential in (P5|), the asymptotic DL length ^ depends, in gen- 
eral, also on k and rj. Starting from the integrable case t] = rjo = 0, and increasing the 
level of classical chaos by increasing 77 at fixed k ^ 1, we find that the sensitivity of ^ to 
an increase in the analyticity of the potential (exhibited by the hopping constants in the 
tight-binding model) is gradually reduced. This phenomenon, which has been considered in 
detail in Sec. VI and is clearly shown in Fig. 7, is a vivid manifestation of classical chaos in 
quantum dynamics. As soon as one reaches a semiclassical regime [r] < 1) of global chaotic 
diffusion {K = Airrik > 1), ^ becomes totally unaffected by any increase in the analyticity of 
the potential. In this regime, ^ is completely determined by the classical chaotic-diffusion 
coefficient D,^^ D/2. 



While QAR-localization is basically different from DL, these localizations are expected 
to look quite similar if 77 is sufficiently irrational and close to r]Q. The numerical results 
shown in Figs. 2-7 strongly support this expectation. In particular. Figs. 5-7 show that in 
cases where the QAR-localization length ^0 is finite (^0 7^ 0) the DL length ^ is well approx- 
imated by ^0 if V is sufficiently close to tiq. In the case of the potential (|5T|), corresponding 
to a pseudorandom Lloyd model, ^ appears to be independent of r], so that ^ = ^o- Using 
then the fact that the KR with the potential (0) exhibits QAR, and ^0 = D/2 exactly, the 
relation ^ = D/2 for this system follows immediately. Ref. ||^ presents a much lengthier 
derivation of this relation, based on the assumption that the pseudorandom disorder can be 
replaced by a truly random one. 



In summary, the presence of QAR in nonintegrable systems is quite useful for study- 
ing several interesting aspects of DL, in particular the transition from DL in local-chaos or 
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almost-integrability regimes to DL in the global-chaos regime, and for obtaining exact lower 
bounds ^0 to the asymptotic DL length ^. 
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APPENDIX A: GENERAL QAR OF PERIOD ONE AND PERIOD TWO 

IN NONINTEGRABLE MKRS 



In Sec. Ill, we have restricted our attention, for simplicity, to a case of period one QAR 
occurring when tj in (13) is a multiple of 27r, tj = 271171 j. Here we shall consider the more 



general case of tj equal to a multiple of vr, tj = mjir, giving both period-one and period-two 
QARs. We believe, but are unable presently to give an exact proof, that this is actually the 
most general case of QAR in the nonintegrable MKRs (|12|). We derive the QAR effective 
Hamiltonian [analogous to ([19|) ] in this general case, and obtain an interesting relation [Rel. 
( |A14|) below] between the values of k^jf associated with period-one and period- two QARs 
in the M = 3 MKR. 

By repeated application of Rel. (J^), we find that for tj = mjir the evolution operator 
(|1^) can be expressed as follows 



U = exp 



M-l 

—ik ^ CjV{9 + rhj7i) 

j=0 



exp(— imovrn^) , (Al) 



where rfij = Y^'^i~^ mg. Consider first period-one QAR. It is clear from ( |Al| ) that this QAR 
may be possible only if mo is even. In this case, one has U = exp[—ikV{6)], where V{6) is 
the function having Fourier coefficients 

Af-l 

Vn = VnY. c.e'"'^^ . (A2) 

3=0 

Here f„ are the Fourier coefficients of V^(^). We denote cj by Cej or Coj depending on whether 
the corresponding rfij is even or odd, respectively. The sum over all Cej (or all Coj) will be 
denoted hj Ce = J2 Ce,j (or Co = J2 Co,j)- Rel. (^) can then be written as follows 

Vn = Vn[Ce + i-lTCo] . (A3) 
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Now, period-one QAR, i.e., U = 1 [without loss of generality, we assume that /3 = in (Q)], 
implies that w„ = for all n. From Rel. ( [A^ ) we then obtain the following conditions for 
period-one QAR, corresponding to two different cases: 

Ce-Co = ifV{e + 7i) = -v{e), 

(A4) 

Ce = Co = otherwise. 
Notice that the second condition [for general V{6)] leads to the trivial result cj = for all 
j if each of the sets {cej} and {coj} contains only one element. This may happen only if 
M = 2 and tuq and mi are both odd. Conditions ( |A^ ) for tj = rrijTc are a generalization of 
condition (p^ in Sec. III. 

If mo is odd, we show that period-two QAR takes place. We find in this case that 
U'^ = exp[—ikV (6)], where V {6) has Fourier coefficients 

2M-1 
3=0 

Here fn'- = J2l^~^ iris, and Cj and rrij are "extended" beyond j = M — 1 by defining, 
for j > M, Cj = Cj-M and rrij = mj_M- From the definitions of fh'- and m^, we see that 
fn'j = mj_M — rfiQ for j > M. Using the last relation in ( |A5| ) , we find that v^^^ = 2t>„ Z^jlo^ '^j 
for n even and v^^^ = forn odd. The conditions for period-two QAR, i.e., [/^ = 1 or v^^^ = 
for all n, are therefore 

Cj- arbitrary if ¥{6 + n) = -¥{6), ^^^^ 

Sjlo^ '^i — otherwise. 
Thus, if V{6 + tt) = —V{6) and tuq is odd, one has period-two QAR for arbitrary values of cj. 
This is a considerable generalization of the period-two QAR in the KR (see Sec. II), discov- 
ered by Izrailev and Shepelyansky [^. The second condition in ( |A6| ) [for ¥{6 + 71) ^ —V{6)] 
is precisely condition (0) in Sec. III. 

We now consider small perturbations of tj near their QAR values m^vr. For definiteness, 
we shall work out in detail here the case of period-one QAR, i.e., tuq even, but we shall show 
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at the end how to extend the results to period-two QAR. Writing tj = mjTr + ej in ([T^), and 
formally expanding in powers of ej as in Sec. Ill, we find, to first order in ej, 

M-l 

^ ~ 1 - E {in^ - k[2iv;i9)h + v;'ie)] + tk^v;\e)} , (a?) 

i=o 

where Vj{9) = J2i=o(^sV{9 + rfisn). Using (p!7|), the expression in ( [A 71 ) can be written, to 
first order in e, as exp(— ieGi), where 

Gi = [n- kVM]^ + k^[AV'i9)f . (A8) 

Here 

VM = E W , [^v'{e)r = E -""^vfio) - k^(^) . (A9) 

j=0 ^ j=0 ^ 

The QE problem in the infinitesimal vicinity of QAR is Giip = gip, and after the gauge 
transformation 

^ = exp[-zkVai9)]i: , (AlO) 
it reduces to the ID Schrodinger equation 



+ k'[AV'{e)Yip = gif . (All) 



The analysis of Eq. ( [AllD is completely similar to that of Eq. (|23|) in Sec. Ill, and the 



conclusions are the same: the QAR-localization length in entirely determined from the 
analytical properties of the function [AV"''''^''^^ 



The results ([A8|) and ( [All|) are valid also in the case of period-two QAR (mo odd), but 
all the sums over j [see ( |A7| ) and ( [A9D ] run now from j = to j = 2M — 1 {tj = tj_M for 
J > M), and Vj{e) is defined as Vj{e) = Ei=oCsV{e + m». 



As an interesting example, we consider the M = 3 MKR with V{6 + vr) = —V{6) 
and rrij independent of j = 0, 1, 2. For rrij = 2, one has the simple case of period-one 
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QAR treated in Sec. Ill, provided, of course, cq + ci + C2 = 0. The value of fee// in ([I9D 
can be easily expressed in terms of two independent coefficients, say Cq and Ci, using (pOD 
- t,)/T = 1/3 for allj]: 

hff = k[2{cl + cj + CoCi)/9]i/2 . (A12) 

Suppose now that the period T is halved, T T/2, leaving all other parameters (includ- 
ing Cj) unchanged. Then rrij = 1 for all j, and, since rfiQ = 3 (odd) one now has period- two 
QAR. The "real-time" period remains then T. Notice also that the condition co + ci + C2 = 
is consistent with ([A6|) [cj can be chosen arbitrarily in this case!]. Using V{6 -|- vr) = —V{6) 
in (|A9D [with M ^ 2M and Vj{e) = Ei=oCsVie + m»], we get 



where 



fcf) =fc[5(cg + c? + CoCi)/6]^/2 • (A13) 



By comparing ( |A13D with ( [A12| ), we obtain the simple relation, valid for all values of Cq and 

Cl, 

= (15/4)1/2 _ ^A14) 
We were not able to discover similar simple relations for other MKR systems. 



APPENDIX B: GENERAL QAR IN INTEGRABLE MKR SYSTEMS 



It is instructive to study the QAR phenomenon in a "linear" version of the MKR, defined 
by the general Hamiltonian: 



36 



H = U + kvie)J2c,ATit-t,] 

^ j=o 



M-l 



(Bl) 



We show that the system is equivalent, effectively, to the linear KR (^), which is inte- 
grable and exactly solvable to a large extent p3[. We then derive general conditions for 
QAR of arbitrary period p in (0), and show rigorously the existence of exponential localiza- 
tion in the infinitesimal vicinity of QAR. 



The evolution operator for (Pl|), from t = —0 to t = T — 0, can be expressed as follows 

M-l 

f/ = exp{—iTjh)exp[—iCjkV{9)]=exp{—iTh)exp[—ikV{9)] , (B2) 

j=0 



where tj = (tj+i — tj)T/T, 



M-l 

V{9) = Y: c,V{e + X,)] , (B3) 

j=0 



Xj = J2i=o'^j for j > 1, and xo = 0- By comparing the last expression for U in ( |B2[ ) with 



Eq. (|TU|), we see that the problem has been reduced, essentially, to that of the linear KR ||2B 
with an "effective" potential V{6). Thus, from the discussion at the end of Sec. II, it follows 
that QAR with arbitrary period p occurs precisely at rational values of t] = t/27i = m/p, 
provided Vsp = for all s. Here Vn are the Fourier coefficients of V{6) in ( p3|) , and are given 
by 

M-l 
3=0 

where Vn are the Fourier coefficients of V^(^). It is now clear from ( |B4| ) that, in contrast 
with the linear KR case (see Sec. II), the requirement Vsp = does not necessarily imply 
the vanishing of Vgp- For simplicity, let us assume in what follows that all r^-'s are equal, 
Tj = t/M for all j, so that Xj = j'^/M. Then, with r = 2TTm/p, the condition Vsp = is 
satisfied for all s if 

M-l 

Cs= J2 9 exp (27ri jsm/M) = . (B5) 

j=0 
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Now, if m and M are relatively prime, the sequence Cs in ( |B5| ) is, up to some rearrangement 
for m 7^ 1, the Fourier transform of Cj (j, s = 0, M — 1). Then = necessarily implies 
that Cj = 0, and the QAR is trivial for generic potentials V{d). If, on the other hand, m 
and M have a maximal common factor m > 1, the condition = implies only that 

m—1 

Cj+rM' = , (B6) 

r=0 

for J = 0, M' — 1, where M' = M/m. To summarize, the M' equations (p^ ) are the 
necessary and sufficient conditions for QAR of period p in the linear MKRs with equal r/s 
(tj = 2'Km/pM) and arbitrary potentials V^(^^). These conditions lead to nontrivial results 
only if m > 1. It should be noticed that the conditions ([B6| ) do not depend on p, i.e., they 
are the same as those for the fundamental (p = 1) QAR with r = 27rm. 



We now consider small perturbations of r near the conditions (P6|). Since the evolution 
operator ( P2| ) is the same as that of a linear KR with effective potential V{9), an exact 
expression of the QE states for irrational r] can be immediately written using the results of 
Ref. m-. 



UO)=^Mmm , (B7) 

where, for all integers /, (pi{6) =16 + 4>{d), and (l){6) is a periodic function with Fourier 
coefficients 

^""^ (n^O), (B8) 



" 1 - e'""^ 

Vn being the Fourier coefficients of V{6), and 0o is arbitrary. The corresponding quasienergy 
is given by uji = vq + It. Now, in the ordinary case of the linear KR, the result analogous to 
( [B8| ) [with Vn replaced by the given coefficients Vn of V{6)] is clearly not defined for n = sp 
and 1] = m/p, i.e., at quantum resonance. In our case, however, with the coefficients Vn 
given by (P^), it is easy to show that the expression ( |B8| ) for n = sp is well defined in the 
QAR limit t] m/p, the limit being taken on some sequence of irrational ?7's converging to 
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m/p. In fact, writing r = (27rm + e)/p and using the conditions (p6|), we find from 
(with Xj = jr/M) that 

M'-l m-1 



j=0 r=0 

After substituting ( |B^ ) into ( JB^ ) and using ([BBD, we obtain 

7,,, M'-l m-1 

hm^ 0.P = E e^-^^-^-'/M' ^ , . (BIO) 

Rel. ( [B10|) shows that the QE states are well defined in the infinitesimal vicinity of the QAR 
of arbitrary period p. We see from (p8|) and (|B10|) that these QE states are localized in L 
space with the same localization length of the (analytic) potential V^(^). In accordance with 
this, the QE spectrum of in the limit of infinitesimal e is pure-point, uji mod 27c = el. 



These rigorous results, which have no analogue in the ordinary case of the linear KR 
p3| , can be derived by an alternative approach, similar to that used for the nonintegrable 
MKRs in Sec. III. The basic evolution operator f/^ can be expressed in the form (|TT|) with 
V{6) replaced by V{9). Then, for pr = 2iTm + e, one easily finds, using (P6|), that is 
given, to first order in e, by exp(— ieGi), where 

Gi = n-0[,(^), (Bll) 

with 0o(^) = liiiie^o 0(^)- In the case of the fundamental QAR {p = 1), with tj = 2TTmj in 
( P2|) as in Sec. Ill, we find that (f)o{6) = kdV{9). Thus the operator ( P 1 1|) may be considered 
as the linear version of (|1^). Its eigenvalues are simply all the integers /, and its eigenstates 
are given by ( pT]) with (j){6) replaced by (poiO). We then see that in the limit e the 
QE spectrum and eigenstates of the approximate evolution operator = exp(— ieGi) agree 
precisely with the rigorous ones obtained above. This may be evidence that the results 
obtained by the self-consistent approach in Sec. Ill are, in fact, rigorous. 
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APPENDIX C: NEARLY-SINGULAR CASES OF 
MULTI-CHANNEL TIGHT-BINDING MODELS 



Clearly, all the hopping constants Wg^n, s 7^ 0, in the singular case [Eq. (^51)] have the 
same range in n (i.e., the range of Wn)- It is thus interesting to consider nearly-singular 
cases of (^) (M > 2) for which the hopping constants W^o,n within a channel are nonzero, 
while Wg^n for s ^ still have approximately the same range in n. The simplest possible 
case where this may happen is when Cj takes only the values 0, ±1, ±2. In this case, the 
Fourier transform Ws{0) of Ws,n assumes the relatively simple form 

where c^^\s) and c^'^\s) are, respectively, the contributions of the terms with cj = ±1 and 
Cj = ±2 to c(s) in (|38|). Using c(i)(0) + c^^\0) = [Rel. (|16D], we find that 

As an example, consider the case of W{6) = kcos{6) — kq, which corresponds to the poten- 
tial ( P7| ) {A = —2/k) and gives a nearest-neighbor ("Lloyd") model ( P^D for M = 2, i.e., 
Wn = except of W±i = k/2 (such a model is studied in Sec. VI). We see from ([Cll) and 



(p2D that there are two interesting limits. If kq ^ max(K, 1), all the hopping constants 
Ws,n are approximately nearest-neighbor in n, and E = c^^\0)ko. If, on the other hand, 
max(K, kq) <^ 1, Ws^n are approximately nearest-neighbor for s 7^ 0, while lVo,n are ap- 
proximately next-next-nearest-neighbor and E ^ c^^\0)kq. The localization length ^ = I/7, 
however, depends only on W{9), and it is always given exactly by Rel. (^81). 

It should be noticed, however, that the choice W{9) = ncos(9) — kq is not a good 
one for obtaining a nearly nearest-neighbor model in the strong-chaos regime {D ^ 1). 
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This is because Eqs. ( p2D and ( ]28|) imply that |/t/Ko| S> 1 in this regime, so that 
Ty(^) = Kcos(6') — Ko = 1 for some 6, leading to a singularity in Eqs. ( PT|) and (|C2D . Let us 
choose, instead, 

2W{e) 



1 - w^{e) 



- iaxi[kV{e)\ = K cos{e) - Ko , (C3) 



I.e., 



It is clear from (|C^ ) that in the strong-chaos regime |iy(6')| ~ 1 for almost all 6. Eqs. (|Cl| ), 
( P2|) , and (|C3|) imply then that the hopping constants Ws,n are all nearest-neighbor in n to 
high accuracy, i.e., Ws^n ~ except of ^^^,±1 ~ c*^^''(s)k/4. The corresponding tight-binding 
model (p6D may thus be naturally viewed as a multi-channel pseudorandom Lloyd model for 
the strong-chaos regime. 
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FIGURES 

FIG. 1. Chaotic orbits generated by 30,000 iterations of the classical Poincare map (|l3|) for a 
M = 3 Hamiltonian (p^), with V{6) = cos{6), cq = ci = 1, C2 = —2, and tj+i — tj = T/3 for all j: 
(a) Local chaos for k = 0.3. (b) Global chaos for k = 0.5. 

FIG. 2. Expectation value of the kinetic energy, £^01 ^ ^ function of the "real time" t = es, 
for the MKR system described in the caption of Fig. |^ with k = 2 and for several values of 
e = 2ir{r] — r]Q) (r/o is a QAR point). The continuous curve corresponds to e = 0.01, the filled circles 
to e = 0.11, the squares to e = 0.17, the filled diamonds to e = 0.25, and the triangles to e = 0.28. 

FIG. 3. Fourier transform Eq{v) of -Eo(0 (^^^ caption of Fig. ||) for e = 0.1 and several values of 
k (see legend). The symbols at the bottom are the theoretical predictions for the peaks positions. 



based on the eigenvalues of the Mathieu equation (26). The peaks for k = 0.01 and k = 1.0 have 



been rescaled by a factor of 100,000 and 10, respectively, for visibility. 

FIG. 4. Steady-state probability distributions {\'ip{n)\'^) over angular momentum nh for the 
M = 2 MKR (co = -ci = 1 and h = T/2) with V{9) = cos{9) and k = 10. The three curves 
correspond to three irrational values oi t] = r/27r: r = 10~^ (solid line), r = 10^^ (dashed line), 
and T = 10~^ (dotted line). The saturation of (|^(n)p) around 10~^^ is due to numerical noise. 



FIG. 5. Same as in Fig. 4, but for the potential (\2T\) with A = 1, k = 1, kq = 0, and 
k = 2. The straight thick dotted line has slope —2j, where 7 is determined from Eq. (p^). For 
r sufficiently small, this slope is expected to be very close to the asymptotic rate of exponential 
decay of {\ip{n)\'^). We observe that this is the case even for r not very small, r = 10^"*^ (dotted 
line). In fact, this system corresponds precisely to a pseudorandom Lloyd model, whose DL length 
^ seems to be independent of r (see Sec. VI and Fig. 8). 
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FIG. 6. Same as in Fig. 5, but for k = 10. The four curves correspond to four irrational values 
of ry = r/27r: r = 10^'' (solid line), r = 10~^ (dashed line), r = 10^^ (dot-dashed line), and 
r = 10^^ (dotted line). The straight thick dotted line has slope —2j, where 7 is determined from 
Eq. dH). 

FIG. 7. Minimal Lyapunov exponent associated with the A^th-order truncation of 
the tight-binding model (|4^ ) (i.e., r restricted to \r\ < N) with k = 10. This 
model is equivalent to the M = 2 MKR in Fig. 4. The several curves cor- 
respond, in descending order at = 15, to the following values of r in (j^)'- 

T = lo-^ lo^^ lo^^ lo-^, lo-^, 2 • lo-^^ 5 • lo^^ 0.1, 0.125, 0.15, 0.2, 0.175, 1.0, 0.25. 

Notice the almost coincidence of the curves for r = 10~^ and r = 10~^, indicating very close 
proximity to QAR. See details concerning the calculation of "Jn in the text. 

FIG. 8. Minimal Lyapunov exponents 7Ar=i2 and 77v=i as a function of r near QAR for the 
same model as in Fig. 7. Notice that 71 appears to be almost independent of r, assuming the 
value corresponding to QAR-localization (r — > 0) for all r. 



47 



PRE - Dana et al 

Figure 2 




PRE - Dana et al 

Figure 3 



0.1 



0.0 



1 








1 


• 
t 
• 
• 
• 
• 
• 

: 

• 
fl 


k=0.01 

k=1.0 

k=2.0 

Ok=0.01;Exact 
■ k=1.0; Exact 
♦ k=2.0; Exact 

L 




♦ o ♦ ♦ ♦ ♦ 
■ ■ ■ ■ 

1,1,1 



0.0 2.0 4.0 6.0 



PRE - Dana et al 

Figure 4 




64 128 192 256 



n 



PRE - Dana et al 

Figure 5 




64 128 192 256 



n 



PRE - Dana et al 



Figure 6 




64 128 192 256 



n 



PRE - Dana et al 



Figure 8 



0.4 



0.2 



0.0 



10 




